Non-exponential relaxation and quantum tunnel splitting in molecular magnet Fe 8 
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Magnetic relaxation in molecular magnets under a sweeping field is studied by taking into account 
local stray fields. It is found that the randomness of local stray field leads to a distribution of 
the relaxation rate which subsequently makes the relaxation deviate from the exponential law as 
predicted by the Landau-Zener model such that the Landau-Zener method needs to be revised to 
deduce an exact tunneling splitting. The tunneling splitting and distribution width of local stray 
fields are derived from the experimental data for molecular magnets Fes- 
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Magnetic relaxation by quantum tunneling in high spin 
molecular magnets has become an attractive field of re- 
search in recent years[l-15]. One of the well-studied sys- 
tems is octanuclear iron(III) oxo-hydroxo clusters Fes, 
which has a well-defined temperature independent re- 
gion below 0.36K.|J g, [§ |, § Theoretically the key to 
understand magnetic relaxation by the quantum tunnel- 
ing is the tunnel splitting. The tunnel splitting and the 
magnetization relaxation are connected via the Landau- 
Zener model. (9[ |l6[ [lTj Since the Landau-Zener transi- 
tion rate is explicitly related to the tunnel splitting, the 
measurement of the change of the magnetization after 
one sweep over the resonant point can give the tun- 
nel splitting 0. Usually the tunnel splitting is very 
small (like the ground state tunneling of Fes), and the 
transition rate due to tunneling is also very low. One 
sweep cannot lead to an observable magnetization relax- 
ation. In this case, the tunnel splitting is determined 
by the Landau-Zener method: || ^, ^| multi-sweeps are 
done and the tunnel splitting is deduced from the ex- 
perimental data in a short time region. Up to now, 
the Landau-Zener method ||, [|, || has served as a ba- 
sic tool to study quantum tunneling in molecular mag- 
nets and many other interesting phenomena including 
the oscillation of tunnel splitting with respect to the 
field along the hard axis and the parity effect for odd 
and even resonance, g H|, IH The tunnel splitting from 
the Landau-Zener method is found to be sweeping-rate- 
independent and agrees with the result found by using a 
square-root decay method.]^, ^, ^ However, there is still 
a puzzle that magnetization re axation under a sweeping 
field shows a clear deviation from the exponential behav- 
ior as predicted by the Landau-Zener model. j| [| Such 
a consequence may lead to a question how the tunnel 
splitting is deduced from the experimental data in short 
time region exactly. This is the main motivation of the 
present paper. 

In the present paper, we start with the biaxial spin 
model with a local stray field to study the magnetization 
relaxation behaviors. It is shown that the "uncompen- 
sated" transverse component of the local stray field leads 
to a distribution of relaxation rates which makes the re- 
laxation in molecular magnets follow a different mecha- 



nism as in some complex systems like spin-glass. j20|, |2lJ 
Our analysis shows that the relaxation is determined by 
two independent quantities: the tunnel splitting and the 
distribution width of the local stray field. Although the 
magnetization relaxation deviates the exponential law we 
can derive the two quantities from experimental data of 
molecular magnets Fes successfully. 

The biaxial spin model for the molecular magnets Fes 
with a local stray field is written asp~8L [191 p2[ 



H = KxS 2 z + K 2 Sl - gu B S (B + h 



(1) 



where K\ > K 2 > , and B is the applied magnetic field, 
h is the local stray field which may originate from the in- 
teractions between the giant spin and the environmental 
spins (including other giant spins or nuclear spins). To 
simplify the problem we assume that h has a Gaussian 
distribution with an equal distribution width in all direc- 
tions 
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Since our main interest is the magnetization relaxation 
under a sweeping field, the external magnetic field is 
taken to be B = {B x ,0,0}: B x = nAB ± ct where n 
is integer, AB is the field interval between neighboring 
resonant tunneling and c = dB x /dt. In the following cal- 
culation, we take K x = 0.321K, K 2 = 0.229K for the 
molecular magnets Fe 8 . || When the field along the easy 
axis is sweeping over the resonant field nAB, the biased 
local stray field h x will be compensated by the sweeping 
field which brings spins into the resonant tunneling and 
leads to a continuous relaxation. If we omit the trans- 
verse component of the local stray field, then all the spins 
will have the same tunneling rate inside the resonant win- 
dow. The resulted relaxation is the simple exponential 
decay according to the Landau-Zener modelpl, p[ 



M(t) = M e 



-n 



(3) 



where T — kP^zc/A and Plz is the Landau-Zener tran- 
sition rate, Plz = 1 — cxp A\j v n c) , A n is the tunnel 
splitting, v n — 2g//sfi(2s — n), k = 2 for n — and k = 1 
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for n — 1, 2, 3, • • •, and A is the amplitude of the AC field 
used in the experiment. In the low transition rate limit, 
i.e., Plz -C 1, which holds for Fes systemj|, the above 
equation leads to 
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In the Taylor series expansion, one has 
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In a short time region, we can keep the first term on 
the right side of Eq. (||) , and then the tunnel splitting A„ 
can be deduced from the magnetization relaxation of the 
molecular magnets. 

However, the tunnel splitting strongly depends on the 
strength of the transverse local stray field. In other 
words, the local stray field will leads to a distribution 
of the tunneling splitting, and furthermore a distribu- 
tion of the relaxation rates. Consequently, the resulted 
relaxation is modified as 



M(t) 



M J dhP(h) e - r(h)t , 



where 



T(h) - [1 - exp (-TrA 2 n (h)/v n c)]kc/A. 



(5) 



(6) 



The tunnel splitting A? (h) can be calculated by the in- 
stanton method, |l9|, |2"3fl 



A2(h)^ 7± (h)A2 Q 



(7) 



where the renormalized factor is caused by the local stray 
field, 7±(h) = [cosh(2qh y ) ± cos(2eJ n /i z )]/2 and ± stands 
for the even and odd resonant tunneling. 
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q = g^ B ^/2{K 1 K 2 
the transverse field, 
we have 
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and A„o is independent of 
In the low transition rate limit, [EL |5l 



M(t) 



Mo J dhP(h)exp{- 7± (h)IV} 



(9) 



where T = fc 7 rA 2 l0 /(^„yl). 

The above equation shows that, in the low transi- 
tion rate limit, the magnetic relaxation is sweeping-rate- 
independent, provided that the sweeping rate is large 
enough, namely, for dB x /dt > 1.0 mT/s.@, || In the ab- 
sence of the transverse local stray field, i.e., P(h) — S(h), 
we have M(t) ~ MqB~ °* for an even resonant tunnel- 
ing and M(t) = Mq for an odd resonant tunneling as 
expected. The effect of the transverse local stray field 
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FIG. 1: Short time relaxation behavior for both odd and even 
resonant tunneling with a = 0.08T and ho = a/A. 



can be observed by doing the integration in Eq.(||). Nu- 
merical results 1 24 as shown in Fig.l show that the mag- 
netization relaxation deviates apparently from the expo- 
nential decay as the distribution width becomes larger 
than 0.03T, which is in a qualitative agreement with 
the experimental measurement. The non-exponential de- 
cay indicates the Landau-Zener model cannot be applied 
to measure the tunneling splitting explicitly. In fact, 
for a system with a distribution of relaxation time, the 
resulted relaxation can be a large variety of shapes of 
decay. |l| The above analysis shows that we cannot 
deduce the tunnel splitting according to Eq. (||) when the 
relaxation deviates strongly from the exponential law. 
On the other hand, consider the local stray fields one 
obtains 
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where (X) = f dhP(h)X. Comparing with Eq.(EJ), one 
can see that the tunnel splitting determined from \M{t) — 
Mo] /Mo under such an approximation is \J (A^) instead 
of A n . From Eq.(Q), one has 



<3±A„o, 



(11) 



where the averaging renormalization factor due to the 
local stray fields is 
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cos(d„h ) 
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(12) 

Dependence of Q± on the distribution width a is plotted 
in Fig. 2. It is seen that Q± increase rapidly with increas- 
ing a and become coincid ent as a > 0.08T. For an even 
resonance, Q + ^ 1, and \J (A 2 ) is always larger than the 
exact tunnel splitting A„o, while for an odd resonance, 
the tunnel splitting */ (A^) is not always quenched once 
the local stray field appears. 
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FIG. 2: Illustration of the dependence of Q± on the distribu- 
tion width a. The dashed and solid lines represent for Q+ a 
Q_, respectively. 



FIG. 3: Relaxation curve of ground state tunneling in Fes 
molecule using a sweeping field with A = 7.2495 x 10~ 2 T 
and dB x /dt = 1.4mT/s. 



Now we are ready to derive the tunneling splitting from 
the experimental data of Fes- The experimental data 
of the ground state tunneling in Fes were provided by 
Wernsdorfer.[[I[ B| To fit the experimental data we choose 
r = 6.5 x 10 /sec and a — 0.05T. Both theoretical 
and experimental results are plotted in Fig. 3. One can 
find that our analytic result fits the experimental curve 
quite well. The estimated value a — 0.05T is consistent 
with the linewidth of the resonance for Feg.Q However, 
the hole-digging method^, |) gives a value 0.03T, which 
is smaller than what we estimated. The tunnel splitting 
from the chosen To is Ao = 7.85 x 10 -8 K, and the renor- 
malization factor is Q+{u = 0.05T) ~ 1.292. Therefore 
we have a/ (Aq) ~ 1.01 x 10~ 7 K, which is very closed 
to 1.0 x 10~ 7 K in the conventional Landau-Zener model 
by Wernsdorfer et al. j|, ^ So what measured in the 
Landau-Zener model isy (Aq) instead of Ao in the short 
time limit. Clearly, Ao and \J (Aq) are two different con- 
cepts. For Fes, since the averaging renormalization factor 
is as large as 1.292, we should consider the effect from the 
local stray field. Our result shows that the magnetiza- 
tion relaxation is determined by the tunnel splitting and 
the distribution width of the local stray field. The later 
quantity leads to the relaxation deviating from the expo- 
nential decay. On the other hand, after the local stray 
field is introduced, the main modification to the kinetic 
equation of the relaxation is to replace A„o with y/ (A^) 
in the low transition rate limit. 

It should be noted that the tunnel splitting from the 
experimental data is sweeping-rate independent since the 
magnetic relaxation is independent of the sweeping rate 
when c > l.OmT/sec. From Eq. (||) we found that, rig- 
orously speaking, the magnetization relaxation depends 
on the sweeping rate c via the expression for T(h). How- 
ever, the larger c leads to a smaller transition rate. When 
the low transition rate approximation becomes valid for 
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FIG. 4: Distribution of ground state tunnel splitting for a Fes 
model with higher order term C(S+ -+- St.) in Hamiltonian 
under a local stray field with a Gaussian distribution. (C = 
-2.95 x 1(T 5 K) 



a sufficiently large c, the transition rate is independent of 
c. This is consistent with the experimental measurement. 
Another problem is that the above analysis is based on 
the adiabatic approximation. It is assumed that all the 
spins have sufficient time to tunnel no matter how small 
the tunnel splitting is. When the distribution width of 
the local stray field is taken into account, the distribution 
of tunnel splitting becomes quite large. Approximately 
there will be a cut-off tunnel splitting A„ c ~ Vncj-n such 
that the spins with the tunnel splitting lower than A„ c 
will have no sufficient time to tunnel. [ fl5| p5| To find out 
the distribution of the tunnel splitting due to the lo- 
cal stray field, we have made Monte-Carlo simulation for 
Fes systems with higher order terms in the Hamiltonian. 
The resulted distribution of tunnel splitting is plotted in 
Fig. 4 for a = 0.05T. The tunnel splitting spreads over 
about 2 to 3 order which is much narrower than that in 
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Mni2 molecules due to dislocation. |l5[ As the local 
stray field from dipolar-dipolar or hyperfine interaction 
in Mni2 is stronger than that in Fes molecules and dis- 
location will lead to a wider distribution of the tunnel 
splitting. The Mni2 system will have a much wider dis- 
tribution of relaxation time than Feg system. It implies 
that the magnetization relaxation in Mni2 system will 
deviate far away from the exponential decay. 

In conclusion we have studied the effect of the local 
stray field on the magnetic relaxation under a sweeping 
field in Fes molecules. The uncompensated transverse 
local stray field leads to a distribution of the transition 
rate such that the relaxation deviates from the exponen- 
tial law. The interplay of the quantum tunneling and the 



distribution of the local stray field determines the mag- 
netization relaxation. Based on this picture we proposed 
a revised scheme to deduce the tunnel splitting A n from 
experimental measurement instead of the conventional 
Landau Zener method. Our conclusion can be general- 
ized to other molecular magnets since the local stray field 
due to the dipolar-dipolar or hyperfine interaction exists 
extensively. 
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